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Abstract 

In this paper we study the Fokker-Plank equation arising in a model which describes the charge and dis¬ 
charge process of ion-lithium batteries. In particular we focus our attention on slow reaction regimes with 
non-negligible entropic effects, which triggers the mass-splitting transition. At first we prove that the prob¬ 
lem is globally well-posed. After that we prove a stability result under some hypothesis of improved regularity 
and a uniqueness result for the stability under some additional condition of the dynamical constraint driving 
the system. 
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1 The model. 

The model introduced in 0), to describe the charging and discharging of lithium-ion batteries, governs the 
evolution of a statistical ensemble of identical particles and is given by the non-local Fokker-Planck equation 

Td t p{x, t) - d x [v 2 d x p(x, t ) + [H'(x) - a (t)) p{x, f)]. (FP1) 

Here H is the free energy of a free particle with thermodynamic state x e R. The probability density p(-, t) 
describes the state of the whole system at the time t, and a reflects that the system is subjected to some external 
forcing. Moreover, r > 0 is the typical relaxation time of a single particle and v > 0 accounts for entropic effect 
(stochastic fluctuation). 

The model 1FP1I has two crucial features that cause highly non-trivial dynamics. First, the free energy H is a 
double-well potential, hence there exist two different stable equilibria for each particle. Second, the system is 
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not driven directly but via a time-dependent control parameter, in our case the parameter is the first moment 
of p, that means we impose the following dynamical constraint 



t)dx = £{t), 


(FP2) 


where d is some given function in time, and a direct calculation shows that IFP2I is equivalent to 

<r(f) = f H'(x)p(x, t)dx + r£(t). (FP3) 

Jr 

The different relation between v and r may cause very different dynamical regimes, which have been studied 
in (5) . We are going to focus to what are called slow reaction regimes, in which we have the coupling 

a 

T — -, (1.1) 

logl/v 


for some parameter a e (0, a cr it)- 

It has been seen in (5) that under the assumption 11.11 . the solutions of IFPli . IFP2I can be approximated, in 
the limit r —► 0 by means of some simpler problems. In particular, during a suitable range of times, the solu¬ 
tions of IFPli . IFP2I can be approximated by means of the solutions of the problem IMS1I - IMS3I described 
later. It turns out that during most of the times the function p can be approximated by the sum of two Dirac 
masses. However, during the range of times in which the approximation of IFP11 , IFP2I is valid, the mass of 
p is distributed in a region with size x of order one. During those times the mass of p is redistributed and, in 
particular, the mass which is initially localized near the point xq is transported to two the neighborhood of two 
different points, denoted as x_,x+. This redistribution of the mass is described by the model IMS 11 - IMS31 
and this is the issue considered in this paper. More details concerning the relationship between the problems 
IFP1I , IFP21 and IMS11 - IMS31 are given in (5j . 

The paper is divided as follows: 

• In Section[2]we give some simple assumption on the potential H appearing in the equation IFP1I and we 
will introduce the problem that we study all along this paper. Moreover we give, in Definition ^. 21 some 
assumptions that the potential H has to satisfy in order that the problem makes sense and is well-posed. 
After that there is some technical result (namely Lemma [231 which determinate some class of potentials 
which are admissible and compatible with the assumptions in Definition l2.2l 


• In Section[3]we prove that there is a unique (in a suitable space) solution of IMS1I - IMS3I in the interval 
(-oo, tg) with to > -oo as explained heuristically in the paper (5]. 


• In Section[4]we prove that the problem is globally well posed in all IK, extending the local result proven in 
Section[3]performed in some interval of the form [-oo, ?o] to the whole real line. 


• In Section^ the last one, we finally prove that, up to sub-sequences (not relabeled) of re-scaled times 
[t m ) m such that t m —>• oo the problem converges to some equilibrium. At the end of such section, namely 
in Suhsection l5.ll we prove as well that such equilibrium is unique and independent by the choice of the 
diverging sequence (f m ) m as long as the dynamical constraint (!{t) satisfies some condition in a vicinity 
of the critical time tz in which the mass splitting transition occurs. 


2 Introduction to the problem. 

2.1 Assumptions on the potential. 

In the paper we are going to assume the following hypothesis on the potential H 
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(Al) H is sufficiently smooth at least Cf (R), and 

H'{x ) = ax+ g(x), 

II^ ,, IIi oo (R) < c< oo, 

with g £ L°° and a > 0. 

(A2) There exist constants x** < x* < 0 < x* < x** and a k < 0 < a* such that 

1. H'(x*) = H'(x**) = a*. 

2 . H'(x*) = H'(x**) = a*. 

3. for each x £ (x*,x*) we have that H'[x) £ (cj*,cr*), H"[x) < 0 and H"{x) ^ 0 for x£ (x*,x*) e . 
In particular the inverse of H' has three strictly monotone branches 


X- 

: (-oo,(T*] 

—■ (-OO.X*] 


: [o-*,cr*] 

—[x*,x*], 

x + 

: [cr*,oo) 

— [(T*,00). 


In what follows we refer to (x*,x*) c as the stable interval, whereas the spinodal region (x*,x*) is called the 
unstable interval. This nomenclature is motivated by the different properties of transport term in IFPlt . In 
both stable intervals adjacent characteristics approach each other exponentially fast, hence there is a strong 
tendency to concentrate mass into narrow peaks. In the unstable interval, however, the separation of adjacent 
characteristics de localizes at an exponential rate in time any peak with positive width. 


2.2 The mass splitting problem. 


Remark 2.1. This is a small remark about notations. 

Suppose are given two functions 

TG:l/cR"-R. 


We write 


if 


F{x\,...,x n ) < G{xi,...,x n ), 
F{x i,...,x n ) < C-G(x i. Xn ), 


for some constant C which is independent from the variables (xi,..., x n ) £ U. 


♦ 


As explained in (5) , at the critical time t2 ~ h we expect that the system undergoes a rapid transition from 
the unstable-stable configuration to a new stable-stable configuration. In order to describe this transition, 
in particular, to predict the mass distribution between the emerging stable peaks, we are going to study the 
mass-splitting model, which describes the peak widening effect (for a full description about the peak widening 
model we refer to (5]) on the rescaled time scale s -[ t - )/t in the limit v — »■ 0. To make the notation simpler 

we will denote the re-scaled time 5 simply t. So, all in all the mass splitting problem consists of the following 
equations 

% = (MSI) 

P [X ’ t] ~ eX P (-pSiy + (1 - m)8 ^- (MS2) 

The asymptotics in IMS2i takes place as long as t-> -oo. The meaning of this asymptotic formula will be 
explained later, in Dehnition l2.2l 

IMS2I codify the necessity to impose asymptotic initial conditions at t — -oo, which have to reflect the fact that 
the mass-splitting process starts in an unstable-stable configuration of a two-peaks model, and that the peak 
on the left is a rescaled Gaussian due to entropic randomness, which means as r —► -oo, se{+,-}, although in 
the following we are going to assume 5 = +, the case 5 = - is simply symmetric, and 



(MS3) 
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this because the transition happens in a timespan so short that allows us to consider the dynamical constraints 
£[t) « £* during the whole transition process. 

The mass-splitting problem hence consists of equations IMS1I . IMS2I and IMS3I . 

Moreover, as well as before we require to have well prepared initial data, i.e. 


a=H'\x o)<0, (2.2) 

b=H"(x + )> 0, (2.3) 

H'(x 0 ) = #'(*+) =<h>e (a*, O, (2.4) 

me (0,1], (2.5) 

mxo + (1 - m)x+ = £*, (2.6) 


where indeed 12.21 - 12.41 codify the unstable-stable configuration in which the process starts, and 12.61 is a 
compatibility condition of the initial data with the dynamical constraint. 


In the following it will turn out to be more convenient to use the distribution of p instead 
itself, to this end we define: 


R(x,t)= f p(y, f) dy. 

J— OO 


of the density 
(2.7) 


Combining IMS It and 12.71 we obtain the following: 


dR 

dt 


[H\x) — a) 


dR 

dx 


( 2 . 8 ) 


Moreover notice that multiplying IMS II by x and integrating in the real line, we obtain, after some integra¬ 
tion by parts and using IMS3I and 12.71 : 

£ , dR 

<j[t)= / H\x) — dx. (2.9) 

J dx 

Integrating IMS2I in the interval (-oo, x) we obtain the formal asymptotic: 

R{x, t) ~ mQ (^r) + (1 - m) X{x>x+} (x), (2.10) 

as t —*• -oo. This convergence takes place in L 1 , although it is uniform in compact sets of R\{x+}. In particular, 
considering a (small) neibourhood U Xo of xo, R (xo + ye a t, t) —■ mQ (y) as t —* -oo in Cf* {U Xo )■ The function Q 
is given by 

1 ff ! 

Q(y)= — e 71 dp. 

V 71 J —oo 

The aim of this paper is to show, under which conditions the transition of the system gives, at the rescaled 
time t = +oo, a new stable-stable configuration, i.e. that there exist two non-spinodal states x_,x+ and a 
numerical value m e [-1,1] such that 


H'(x-) = H\x+), 


(m_ - m)x- + (m+ + m)x+ - £* - m_x_ + m + x+. 


The existence and the uniqueness of m is not obvious since the mass-splitting problem involves two sub¬ 
tle limits. First one has to show that the asymptotic condition IMS2I gives rise to a well-posed initial value 
problem at t = -oo. Second, one has to guarantee that solutions do not drift as t —> oo along the connected 
one-parameter family of equilibrium solutions. 

We need some general assumption on the function H in order to prove that the problem is well-posed. 


Definition 2.2. We say that a function H e C? (IR) satisfies the condition H if, for any function <p(t) = cr(t)-cro e 
C°°(-oo, T ), T e R satisfying |</>(f)| Me _(2a+5K ,<5 > 0 and any K e IR, there exist a unique solution to the ODE 


problem 


( 2 . 11 ) 


X(t,K) = x 0 + Ke- at +Y(K t f), as - oo. 


4 



















With Y{K, t) = o(e (_a+5)f ). 

Moreover for every fixed f e (- 00 , T] the transformation 

KeM~ X{t,K), 


is a one-to-one transformation of the real line in a interval (X_(f),X|.(f)) where the functions X ± ( t) solves the 
problem 

< ( 2 . 12 ) 
X+(f) — x+ + Y+(f), as t —- 00 , 

With Y+(f) = o(e (_a+5)t ). 

Where H'{x 0 ) = H'(x±) — (Tq. Finally for any fixed t e (- 00 , T] we have 

lim X{t,K) = X ± {t). 

K—*± 00 

Definition 12.21 is vacuous if the set of potentials satisfying such condition is empty. For this reason the 
following lemma gives an example of a class of potentials that satisfy the Condition H. 

Lemma 2.3. Everyfunction H € C 2 o (R) such that satisfies the compatibility condition H' (x_) = H'(x 0 ) = H ! (x + ) = 
do satisfy the conditions in Definition \2.2\ 

Proof. We are going to subdivide the proof in several steps 

Step 1 We want to show that the solutions of 12.121 are well defined. We are going to do this for the characteristic 
X+, the other case is similar. Let us write X+ in the following form 


X + (t) = x + +Y+(f), 

where Y+ is considered to be a perturbation. Than the equation 12.121 reads as 

Y; = -bY + + (t>(t)-0{YZ), 

and hence if such a solution exists than it has to take the following form 

r+(f)= f e“ 6(f “ s) [0W-O(y + 2 W)]ds. 

J— OO 

Define the following operator 

r + [F + ](t)= T e- bit - s) l(P(s)-0{Y^s))}ds, 

J — OO 

we show that T + admits a fixed point in the space 

/(t 0 ,<5) = {/:|/(t)Ke ( - n+5)f ,t< t 0 ,8 > oj, 


endowed with the norm 


UWnt 0 ,8) = sup-!|/(r)|-e (fl 5)f }, 

t<t a 1 > 


(2.13) 


^ ae~ [2a+m + c 2 e 2{ ~ a+5)t < e ( - n+5 h ; 


for to sufficiently negative and 8 small. 

To do so take Y+ e J(to,8) and evaluate 

|r + [Y+](t)| = I r e- b(f - s) [</)(5)-0(F + 2 (s))]ds 

I J —OO 

for t sufficiently negative, since c,- = c; (M, a, b, 8 ). 

At this point take Y+ ,i, Y +j 2 e /(to,5), and consider 

|r + [F + ,i](t)-r + [F+, 2 ](t)| < T e- Ht - s) \0[Y 2 tl {s))-0[Y 2 2 ^)\ds. (2.14) 

J —00 
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In particular we point out that the functions 0\Y+. (s) J have an obvious explicit expression, as explained 
in equation 13.71 , in particular, in our case 

0(y + 2 ,.(5)) = p +ii (f) = H’(x + + Y +it m - [H'(x + ] + bY +/l (t)). 

With this consideration we can rewrite equation 12.141 as 


|r + [r + ,i](f)- r + [Y+, 2 ](f)| < 

[‘ e~ Ht ~ s) (| H\x + + F + ,i(s)) - H\x + + F + , 2 (s))| + b|F + ,i(s) - F + , 2 (s)|) ds. (2.15) 

J —oo 

In particular, since H e Cj 3 c (M) we can say that H' e Cf' 3 (R), which means that 

| H\x + + F + ,i(s)) - H\x + + F + , 2 (s))| < L | F + ,i(s) - Y + , 2 (s)| , 


for some L > 0 that depends only on H and on a compact Kc IR sufficiently large. 

With this consideration we can rewrite 12.151 as 

r ^ 

| t + [y + ,i] (r) - t + [y +i2 ] (r) | < c VC, h, b) / e~ Ht ~ s) \ y + i1 w - f + , 2 (s) | d s 

J —oo 

< II Fo,i - Vo, 2 II J[ta,8) -e~ bt C e [b ~ a+S)s d 5 < e^ a+5) r ■ II Yb.i - F 0 ,2 II J[to , S ), ‘ 

J—OO 

proving that it is a contraction for t ^ to and concluding the proof of the first step. 

Step 2 We want to prove that for any fixed K there exists a t* = t* (K) such that 12.111 has a unique solution. 

To prove this we write 

X{K, t) = x o + Ke~ at + Y{K, t), 

the proof is the same as in the first step only considering Y ( K, t) belonging to the following space 

Kt 0 ,S,K) = {/ e C°|/(f)| < C K e [ ~ a+S) r , t < to,6 > o}, 

with p defined as in 13.141 . with the norm 

ll/ll j(t 0 ,8,K] = sup||/(r)|-e ta " 6)f l. 

t<t 0 1 1 

We point out the fact that, a priori, could happen that limjf^ ±00 t*{K) = -oo. 

Step 3 In this step we prove that there exist a time ~t — t(<p, <5) such that the solutions of 12.111 are defined uni¬ 
formly for all t ^ ~t. 

To do so fix a K that we might as well consider to be positive. Since as long as \X'(t, K) \ < oo we can extend 
the solutions of 12.111 to the right off, see (2] Theorem 1.1, Chapter 2] for such extensibility theorem, we 
suppose than that there exists af = ~t{ K) such that 

X\t,K)^^oo. (2.16) 

For instance X'{t,K) ~ ^ + 00 . The case X'(t, K) ——- -00 is similar. 

If this happens, than, since X' - cr 0 + (pi t) - H ! (X) —► +00 then it must be true that 

-H'[X'(t,K))l-^+ 00 , 
and hence considering the structure of H' 


X'{t,K) 


t^t 


-00, 


which contradicts the hypothesis 12.161 . 
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Step 4 What is left to prove in this last step is that 

lim X{t,K) = X+(t). 

K —++00 

First we claim that for every t ^ t there exist a K > 0 such that X'{K, t) > 0. As usual for K < 0 the 
procedure is similar. 

Indeed as we have seen in the second step for every rescaled time t there exists a k small enough such 
that the following approximation holds 

X(K,t) = x 0 + Ke~ at +o[e~ at y (2.17) 

But then 

X\k, t) = o{t) - H' (X(k, f)) = -aKe- at +o[e ~ ar ) > 0. 

Where a is defined in 12.21 . Since the approximation in equation 12.171 is valid for every k as t —► -oo 
we can say that limj__ 00 K{t) - +oo. Now, considering the flux structure of 12.111 . we can argue that 
X'{K, t+) > X'{K, t) in some right neighbourhood of t. I want to show that, given an 77 > 0 that we are 
going to consider small, such that X'[K, t) > r), the system reaches in a finite time the state 

X'[K,t) = ri, (2.18) 

where it is not restrictive to set 

t- inf{f+ > t: X'[K, f+) = 77 }. 

We remark the fact that, for every 77 > 0 and small there exists a state satisfying 12.181 , this because if we 
suppose X’ ( K , 0 for every t ^ t we would obtain that 

X[K,i)-X(K,t)>fi[i- t), 

hence, since X + (f) is bounded (see Step 1) there would exist a finite time t such that X ( K, t) = X+ (7), but 
this contradicts the uniqueness of solution for 12 . 111 . 

We want to see that under these assumptions X[K, s) reaches the configuration X {K, t) in a finite time. 
Since the characteristics are uniformly bounded as we have seen in the third step we can say that X ( K , t] - 
X ( K , t) ^ Co, independently of t,t. Moreover 

X[K,t)-X{K,t) = J X’(K,z)dz=[t- t)X’(K, t). 

Where in the last equality we have used the mean value theorem. Moreover 

X'[K,t)^ inf X'(K,() > — >0, 
feffi] N 

since X'(K,s ) > 0 in [t,t\. Whence (t- t) < < oo. Once we have proved that such a state can be 

reached, in a possibly large, but finite time, and since all this procedure has been done independently 
by the starting time f we can state that, taking t sufficiently negative there exists a t such that 

H’ [X [K, t)) = (Jq + (f) (f) - 77 . (2.19) 

But <p (f) is an o f j function, hence we can assert that for t sufficiently negative 

3 n 

00 - -p < ^0 + <P (tj - V < (J 0 - -, 

whence 

<r 0 + 0 (f )-77 = 0-0 + 0 ( 77 ). ( 2 . 20 ) 

At this point, since X ( K , t)> X (K, t), and since the state X ( K, f) has to satisfy 12.191 . considering more¬ 
over 12.201 we can say that X ( K, t) - x + + O ( 77 ), and, thanks to Step 1 we can assert X+ (f) = x + + o(l), 
whence 

|X(r, t) - X + (f)^ 0 ( 77 ) + o(l), (2.21) 

and indeed the right hand side of 12.211 can be made as small as we want, concluding the proof. 

□ 
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3 Well posedness for t ->• -oo. 

We want to see that there effectively exist some p satisfying IMS1 1 , 1MS21 , formalized we want to prove the 
following proposition. 

Proposition 3.1. Suppose H satisfies the conditions in Definition s. 2\ andH'" absolutely continuous in compact 
sets ofR. Suppose also that the following stability condition holds: 

[1- m)a+ mb> 0 (3.1) 

with a and b defined respectively as in 12.21 and 12.31 . For any mE (0,1] thereexistaT £ R and a unique solution 
R of 12. 8 ) , 12.91 . monotonically increasing in x for any t e {-oo,T] such that R(x, t) e [0,1] and satisfying 12.101 , 
i.e. 

\\R(‘,f)-mX{x>xo}-(l-™)X{xZx + }\\ L i m -0, (3.2) 

as t—> -oo. 

Proof. In order to prove the existence of solutions we define the following class of functions: 

K m ,t(S) = { (TeC Q \(T{t)-(Ja\^Me- {2a+5)t , fe(-oo,T] }, 
for some 8 > 0 and sufficiently small. 

Our goal is to show that for T < 0, \T\,M sufficiently large it is possible to define a transformation from 
Km,t(S ) to itself, whose fixed point is equivalent to solving 12.81 . 12.91 and 13.21 . 

Given a e Km,t{8 ) and K e R let X{s,K,cr ) = X[s,K ) and X ± {s,o) - X+(s) respectively the solutions of 12.111 
and 12.121 . which are well defined due to the fact that FI satisfies the conditions in Definition l2.2l We define 
then R (x, s, cr) as follows: 

f R{X{t,K,(j),t) = mQ{K) + [l - m)X{ x ^ x+ (.t,a)} (X[t,K,o r)), 

< R{x, t,o) = 0, if x<X-(t,a), (3.3) 

l R{x, t,o) - 1, if x>X + (t,cr). 

Our goal is to obtain a function cr e Km,t(8) such that satisfies equation 12.91 . with R defined by means 
of 13.31 . We want to obtain a linearised version of 12.91 where it will be implicitly assumed thatjj can be 
approximated for t —>• -oo as two dirac masses at Xo(t,cr) = X(t,0,cr) and X + (t,o) respectively. To this end 
rewrite 12.91 as: 


rX+(t.a) gp rX+(,t,a) gj{ 

a(t) - H'(Xo(t,a)) — [x,t,a) dx+ / \H'(x) - Fl'(Xo(t,cr))] — (x, t,u) dx+ (1 - m)H'[X+{t,o)), 

JX-(t,a ) dx JXAt.a) dx 

whence we obtain 

rX+{t,o) g R 

a(t) = mH\X 0 lt t <j)) + a-m)H\X + (t,(j))+ / \H'[x) - H'[X 0 {t,o))\ —(i,t,ff)dx (3.4) 

Jx-R,o) dx 

At this point, using Taylor expansion we may write 


H'{X 0 {t,a)) = H'{x 0 ) + H"{x 0 ) (X 0 (t,cr) - x 0 ) + p 0 (f,cr), 
H'{X + ( t ,cr)) = H'{x+) + H"{x + ) [X + {t,o)-x + ) + p + (t, cr), 

where indeed po,p+ are reminder with the obvious explicit expression: 


po(bcr) = H'{X 0 [t,a)) - ( H'{xq ) + H"[x 0 ) (X 0 (r,cr) - x 0 )), (3.6) 

p+(t,cr) = H'{X+{t,a)) - ( H'{x +) + H"(x+) (X+(t,cr) - x +)), (3.7) 

hence, considering 13.51 into 13.41 we obtain the following expression for cr(f): 


<j{t)- [mFt'{x o) + (1 - m)H'{x+) \ + ma[Xo(t,cr) - xq ) + (1 - m)b{X+[t,a) - x+) + 

dR 
dx 


rX+(t,a) gp 

/ \Fl'(x) - //'(Xothcr))] — (x, t, cr) dx+ mpo(t,cr) + (1 - m)p+(t,cr), 

JxM.a) 
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with a, b as always defined by means of 1 12.21 , 12.31 . 

Notice at this point that [mH'{x o) + (1 - m)H'{x+)\ = <jq, set <p — cr- erg and set 


](t, a) = 


L 


X+(t,or) gR 

\H'{x)-H'{X a {t,(j))\ —(x, f,cr)dx, 

X-(t,a) dx 


(3.8) 


with this consideration we can write equation 13.41 as 

<p{t) — ma[Xo[t,a ) - xg) + (1 - m)b{X+{t,o ) - x + ) + mpo{t,o ) + (1 - m)p+(t,c r) + /(f,cr). (3.9) 


Our next aim is to express /(f, cr) in a more suitable form, to do so integrate by parts equation 13.81 obtaining 


that 


rx+(t,a) , 

/( t, cr) = - I H" (x) R (x, t, o) dx + R (X+ ( t, a), t, cr) \h' {X + ( t, a) - H 1 {X 0 ( t, a) . 

JX-(t.a) [ J 


)X-(t,o ) 

Moreover, using the following change of variable 

x=X{s,K,a), 

dXCs^o) 

dK 

and considering the Condition (ff), that ensures that the intervals e P is transformed by X(s,-,n) into (X_(s,cr),X + (s, cr)) 
we obtain, using moreover that R (X + ( t, a) , t, cr) = mo 


Kt,(j) = - [ H" {X{t,K,a))R{X{t,K,a),t,a) 

Jr 

and, using also 13.31 

./ 


dX{s,K,a) 

dK 


dK++m 0 [H'{X+(t, a))-H'(X 0 {t, a))]. 


. „ dX{s,K,a ) . , , . 

J{t,o) = -mo j H"(X(t,K,a))Q(K) - dK+ m 0 \ H [X + (t,a)) - H (X 0 {t,a)) \ 


= -mo f 4~{H'{X{t,K,a)))Q{K)dK + mo\H'(X + {t,a))~ H'(Xoit.a))] 
Ju OK L J 

= m 0 f [H' (X(t,K,<r)) - H' (X+(r,cr))] Q'[K)dK 

JU 


Obtaining hence, in the end that 

Kt,cr) = m [ [H\X{t,K,<j))-H\Xo{t,(j)j) Q'(K) 
Jr 


dK. 


(3.10) 


At this point we want to linearise equation 12.111 and 1 12.121 using Taylor expansion for H' around xq and x+, 
obtaining the following equations 

d 

— (X 0 (t,cr)-x 0 ) = -a(X 0 (t,cr) - x 0 ) + cp(t) - p 0 (t,cr), (3.11) 

df 

^-(X+{t,o) - x+) = -b(X + (t,<r) - x+) + (p{t) - p+{t,(j), (3.12) 

df 

d 

— X{t,K,<j) — -a(X[t,K,<j) - xo) + <b{t) - p{t,K,a) } (3.13) 

df 

with p in 1 13.131 reminder of the form 

p(t,K,a ) = [H’{X{t,K,(j))-H'{xo)-a{X{t,K,G)-xo)], (3.14) 

and not a probability density. 

We want to linearise also equation 13.101 , to do so write first 

0 H'{X{t,K,a )) - H'(Xo(t,a))) = [H'{X{t,K,a)) - <j 0 ) - [H'(XoU,a))-a 0 ), 

= [H'{X{t,K,a)) - H'(x o)) - [H'(X 0 (t,a)) - H'(x 0 )), 


9 













whence we obtain 


J[t,u) = ma f (X(t,K,cr)-x 0 )Q'(K)dK-mb [ (X + (t,a)-x + )Q'(K)dK+ 

J re Jr 

+ m[ p{t } K,o)Q\K)dK-m [ p+(t,cr)Q'(K) dK, (a) 

Jr Jr 

and also, considering that / Q'(K ) dK = 1, 

J{t,a)-ma\ {X{t,K,<j) - xq)Q\K) dK-ma{X(,{t,<j) - xq) + m f p(t,K,cr)Q'(K) dK-mpo(t,cr). (b) 

Jr Jr 

We remark the fact that (a) and (bJ are different linearisation of J{t,<j). 


On the other hand we can remove the leading order in the asymptotic of X[t,K,u) - xq in order to rewrite 
the problem as a linearised one plus a perturbative term. We write: 

X{t,K,(j)-x 0 = Ke~ at +Y(t,K,cj) } 

Xo(t,cr) - xq = Y 0 {t,a), (3.15) 

X+{t,(j)- x+ = 7+(f,cr). 

Substituting in 13.91 the results in 13.151 we obtain 

<p{t) = m[H"(x 0 )Y 0 {t,cj) +p Q {t,(j)] + (l-m) [H" (x+)F+(f,cr) + p+{t,a)] + /(f,cr). 

Considering moreover the result in equation © 


(p{t) = m [H"{xq) Yo(f,cr) + p 0 (£,cr)] + (1 - m) [H"{x + ) Y+{t,o) + p+(f,cr)] + 

ma f Y(t,K,cr)Q'(K)dK-maY 0 (t,cr) + m [ p(t,K,a)Q'(K)dK- mp 0 (t,cr). 

t/[R J |R 

So, at the end, considering that / Q’dK - 1 we obtain 


< p{t ) = m[aYo(t,o) + po(hcr)) + (1 - m) [bY+{t,o ) + p+[t,a)) + 

+ ma f (Y(t,K,cr)-Y 0 (t,cr))Q'(K)dK+m f (p(t,K,a) - p 0 (t,a))Q'(K)dK, 

Jr Jr 

and indeed considering 13.151 and linearising using Taylor expansion we obtain: 

d 

— Y 0 (t,cr) = -aY 0 (t,cr) + (p(t)-p 0 (t,cr), 
dt 

d 

— Y+{t,o ) = -bY+{t,o) + (p{t) - p+(r,cr), 
df 

d 

— Y(t,K,a) = -aY(t,K,cj) + <p(t)-p(t,K,o). 
dt 

13.171 - 13.191 are simply 13.111 - 13.131 with the considerations in 13.151 . 

At this point we define the following linear operators: 

nt nt 

U 0 f{t) = / e~ a[t ~ z) f{z)dz , ZJ+/(f)= I e~ b[t ~ z) f{z)dz. 

J—o o J— oo 

Apply the variation of constants method to 1 13.171 — f3TT9l . obtaining 

Yo(f,cr) = U 0 (p{t)-Uap 0 {t,o), 

F+(f,cr) = U+(p{t)-U + p + {t,c r), 

Y (t,K,a) = Uo(p(t)-U Q p{t,K,o). 

Using these equation we can transform 13.161 into: 


(3.16) 


(3.17) 

(3.18) 

(3.19) 


< p{t ) = m[a(U 0 (p{t) -U 0 p 0 {t,(T)) + p 0 (t t cr)] + (l-m) [b{U+(p(t)-U+p+{t,(j)) + p+(t,cr)] 

+ ma\ (Uopo[t,cr)-Uop[t,K,a)) Q'[K)dK+ m f (p(t,K,cr) - po(t,a)] Q'(K) dK. (3.20) 
Jr Jr 
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Rearranging the terms we obtain the following equation: 

C(p{t) = W{t) 


(3.21) 


where 


and 


C t nt 

C(p{t) - cp(t) - ma \ e~ a(t ~ z) cp(z)dz-(1-m)b e~ b(t ~ z ^(p[z) dz, 
J—oo J—oo 


(3.22) 


W(f) = m[-aU 0 lp 0 ]{t,a) + p 0 (t,<j)) + [l-m) [-bU+[p+]{t,(j) + p+(t,o)) + 

+ ma f lUolpo\{t,<j)-U 0 lpKt,K t <j))Q\lC)dK+in [ (p(t,K,cr) - p 0 (t,cr)) Q'(K) dK. (3.23) 
Jr Jr 

As it is shown in Lemma lA. 1 I t he operator C defined in 13.221 can be inverted. Moreover, as proved in Lemma 
IA.2I we can explicitly solve the equation 

£G(x) =8 a {t), 

whose solution is given in IA.31 . 

Thanks to the considerations above, using the Green function G we can write the solutions of 13.211 as 

(pit) - f G(tj- t)W(rjJdT) (3.24) 

J—OO 

we want to solve this equation using a fixed point argument. To do so recall the space 

K M ,to ^ = {</> e C ° : l CT « - CT ol = 10(01 < Me~ {2a+S) t ,MeU,te (-oo, t 0 ]}, 

Endowed with the norm 

11011 K ut iS)= sup m)\-e [2a+S)t , 

M ' r ° t€(-OO,t 0 ] 

we want to apply Banach fixed point theorem for the operator 

T[(j)][t)=f G(rj-t)W{Tj)dr), (3.25) 

J—OO 

where the function W ( t) — i.e. there is a direct non-linear dependence on the function <p. To do so we 

proceed as usual in two steps, 

1. Check that T maps Km,t(8) to itself for to sufficiently negative. 

2. Check that T is a contraction. 

Indeed in order to verify |T] and [2] we need some estimate for \ W(t)\, which are given in Proposition IA.3I i.e. 
equation IA.9I tells us that 

\W(t )\< Ce~ 2at , 

for some positive C-C(a, b, m,8) < oo uniformly in 8 and t < to sufficiently negative. 


This is the first ingredient in order to prove the fixed point for T. 
Since G £ K m j ( 8 ) for each 8 > 0, then 


\Tmm- 

and indeed 


I f 

J —OC 


G(rj- t)W(T))dT) 


^ CMe [2a+S) 




e - (4fl+5) t?d n^—^Le-^t-. 
, 4a+ 8 


CMe‘ 


.St 


4 a + 8 


-(2a+S) t 


Ce‘ 


St 


< M, 


4 a + S 

for t sufficiently negative. This proves that the operator T effectively maps K m j 0 {8 ) onto itself if to is suffi¬ 
ciently negative. 


11 










To prove that T is indeed a contraction on K m j 0 (5) we will have to repeat some calculations which are 
made explicit in Pronosition lA.41 

Considering the bound given in 1A.101 we can evaluate 
|n<W(f)- T[(p 2 ](t) | = I f G(tj- t) [W{q,(pi)-W(q,(p 2 )]dq 

\J—o O 

r t rt 

< |G(f7-O||[W’(i7,0i)-W r (i/,<!&2)]|di/<||0i-02|| jr \G(q-t)le~ (2a+S) Uq 

J-oo M 'h J-oo 

GeK m tor t 

< IWi-Mk M t (S)f e- (4fl+25) ^t 7 <||0 1 -^ 2 || 5) . e - 2 ‘ 2a+5)f . (3.26) 

M ’ r 0 J-oo M - l 0 

Multipling both sides of 13.261 for e (2fl+l5) t , and taking the sup for t f to we obtain 

|| ri^Kf) - T[<p 2 \{t) 1^^ < e- (2fl+5) \\<l>i ~ <p 2 , 

which guarantees that T is a contraction concluding the proof of the theorem. □ 


4 Global well posedness. 


From now on the re-scaled time t is going to be called t. 

The following lemma is the starting point of our analysis. 

Lemma 4.1. Consider a as defined in 12.91 . than, a is completely determined by the evolution of the character¬ 
istics X±. 

Proof. Apply integration by parts obtaining 


/■x+Cfl 

<j(t) = H'{X + {t)) - H"{x)R{x, f)dx. 

JXM) 


In the same way we can obtain the following equation 

r-X+W 


CA+UJ 

o{t) = H'{X-{t))+ / H"{x){l-R{x,t))Ax. 

JX-tt) 


(4.1) 


(4.2) 

□ 


In view of the simple computations in Lemma l4~Tl it is equivalent to prove global existence for a or for the 
characteristics X ± . Recall that 

( x' + m = -H'[x + m+<j(t), 

1 X'_(t) = -H'{XM)) + (J(t), 

considering 14.U and 14.21 we obtain the following new differential equations for the characteristics 


< 


rA+UJ 

+ (£) = -/ H r '(x)R(x,t)dx, 

rX+W 

’_{t) - I H"{x) (1 - R(x, t)) dx. 

JX-(t) 


(4.3) 


Which is very interesting since we see explicitly in 14.31 that the evolution of X+ is not influenced by a, hence 
we can write 14.31 as X’ — f{X, t ) with /: R z x H ^ R 2 and has the following explicit formulation 


( r xi r xi 

fix, t) — [ - I H"{x)R{x, t) dx , 1 H"(x){l-R(x, f))dx 

V JX2 JX2 

At this point we can start to study the system 14.31 . 

Proposition 4.2. Under the assumptions 12.11 the system 14.31 admits unique solution globally in R. 


(4.4) 


12 








Proof. The function / is Lipschitz continuous and sub-linear in the variable x, hence apply (7j Theorem 2.17] 
obtaining our global result. □ 


The global existence statement proved in Proposition l4.2l gives. thanks to the computations performed in 
Lemma fOl a global existence result also for the function cr, which we showed that the evolution of o can be 
completely described in terms of the evolution of the characteristics X+. 

We would like to refine even further this result, namely we would like to be able to give some L°° bound for the 
characteristics X ± and, consequently, for the function o that we might need in the future. 


Lemma 4.3. Define AX{t) — X + (t) - X_(t), and consider a potential H which satisfies 12.11 . In this case AX e 
L°° (R). 

Now we can obtain the L°° estimates for a,X+. 

Proposition 4.4. Consider a potential H satisfying 12.11 . than we have that the function a defined in 14.11 or 
equivalently 14.21 and the functions X+ described by the system 14.31 are not only defined in all [R but they also 
belong to L°° (R). 

Proof The proof is, at this point, very short. Consider a as in 14.11 . namely 

r x + {t) 

a{t) = H\X+W)- H''{x)R{x,t)Ax. 

JX-W 

We know thanks to the previous lemma that |AX| < L hence | ^ H" (x)R{x, t) dx| ^ || //|| • L < oo. Which 

means that f*H"(x)R(x, •) d x e L°°(IR). Considering 12.11 . H'(x) = ax+ g(x) where g is a L°° function. 

With these considerations we obtain that 

<j{t)- J H r (x)pdx-aJ xpdx + J g(x)pdx ad* + llg'llz. 00 (iR) < oo 

We’ve obtained hence that cr,f^ + ^ H"{x)R{x,»)dx e L°°, which means, considering 14.11 . that H'(X + ) e 
L°° which implies that X+ e L°°. 

To prove that X- e L°° the reasoning is the same considering a as in the equation 14.21 . □ 


5 Stability for t —► +oo. 


Proposition 13.11 assures us that as long as -oo the solution to IMS 11 stabilize to a convex combination 
of Dirac-d measures localized in two points. We expect to recover, after the mass splitting process, a similar 
configuration. 

To do so we are going to show that the system effectively converges to such a form via a standard stability 
argument in a neighbourhood of t — + 00 . 


Lemma 5.1. Define the following function 

E t — J H(x)pdx<oo, 

where p — pit, x) solution of 1MSU . Then we have thatE t is decreasing in time. 


(5.1) 


Proof. Derive equation 15.11 , integrate the obtained equation by parts and apply Jensen inequality obtaining 
that 

H(x)pdx\ 0 


dt 

Which is exactly what we wanted to prove. 


/■ 


□ 


At this point it is clear that as t —► +00 we have that the system stabilize to some value E +oa f E t for each 
t e [R u {- 00 }, and that 15.11 holds. 
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Theorem 5.2. There exists some sequence sequence (t m )meN such thatt m m °°» +oo. Setp m {x) = p { t rn , x), than 
there exist some X-, xq, x + , such thatx+ £ [H" ^ 0} andx o £ { H" < 0} and 

PmW — ]T mi6 Xi {x) (5.2) 

ie{—, 0 ,+} 


with ni- + mo + m+ - 1 and 

H\x-) = H'{x + ) = H\x o) 

Proof. To prove the theorem proceed as follows. Define the dissipation of the system as 

-D(f) = —| J H(x)pdxj = - J (H'(x)) 2 pdx+H'(x)pdx j ^0 

From this definition we obtain that / ^ D(t)dt — - -f^E t dt= E-qq-Et < E-„q. 

At this point, hence, we have obtained that 


f 

Jr 


D(t)dt < oo, 


(5.3) 


(5.4) 


(5.5) 


which implies that there exists a sequence (t m ) rn such that D[t m ) —* 0 as m —* +oo. In particular let us select 

m—*0 O , , _v r \ /-v • 

some sequence t m -* +oo, we have that D(t m ) —* 0, i.e. 


J [H’{x)) 2 p m dx-H'(x)p m dx 


0. 


(5.6) 


Where p m {x) =p{t m ,x). 

Moreover the set [pit, •)) is bounded in rca(Vi), the set of regular measures, with finite variation. By Banach- 
Alaoglu theorem is weak-* compact, hence there exists a (not relabelled) sequence of diverging times ( t m ) m 
such that p m p. 

15.61 verifies if and only if H' — k for some keU p-almost everywhere, for p - * - lim m p m . This consideration, 
together with 12.11 . and the fact that H' is strictly monotone in its invertible branches, implies that 

Pm ^ tniSxjt 

/£{-, 0 ,+} 

where the points xt,ie {-,0, +} have to satisfy the condition 15.31 since H' - k p-a.s. or, else p m 8 X for some 
x £ [x*,x*] c . 

A priori could as well happen that p m —- 8^, we want to exclude this eventuality. 

To do so consider cr{t m ) — f E['{x)p m {x) dx —► oo but a e L°° as proved in Pronosition l4.4l hence this is absurd. 

□ 

Proposition 5.3. Set a sequence (t m ) m such that p,„ — Xi€{-,o,+) m iS Xi ■ Than for any f £ C[D) where D is the 
following compact set of U: D - [D-,D+] with 


D- = inf X-it), 

teR 

D + = supX + (f), 

teR 

where X+ are solution of 14.31 . the following equality holds 

lint / f(x)pit m ,x)dx- V mifixi). 

Proof. We point out at first that \D± \ < oo thanks the fact that X± e L°° as shown in Prooosition l4.4l 
With these considerations, for every time t the mass of the entire system is concentrated in the set [ X- ( t) , X + ( t) ] , 
this has been described in detail in Lemma [231 considering also the global result in Prooositior l4.2l This im¬ 
plies that the support of the probability density p is compact, and independent from the time t, namely the 
set D. The claim hence follows since, in this setting, [C° (D))* = rca[D) and pit,-) e rcaiD) for all t £ IR. □ 
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Corollary 5.4. Set a sequence {t m ) m such that p m — o,+) m iS Xi ■ Then 

lim a{t m ) = H'{Xi) i G , 0 , +} 

m—*o o 

e= Y m ‘ x i 

z e {—, 0,+) 

Proof. AddIv ProDosition l5.3l to the functions crit) - f H'(x)p(t, x) d x and - f xp ( t, x) d x obtaining 

<J{t m ) = Y miH\xi) + o{\) = H\xi) + o{\). 
ie{-,0,+} 

and 

( = I xp m {x)dx - Y m i x i> 

J z'e{-,0,+} 


(5.7) 


(5.8) 

□ 


Corollary 5.5. For any teU and any f eC[D), it is true that /e L 1 (IR, p(t, •)). 

Proof, f f{x)p{x, t) dx ^ maXjED |/(x)| < oo. □ 

Lemma 5.6. Select a sequence (t m ) m such thatp m -^L;e{-,o,+} uvS Xl ■ Supposexi x+,x*. Define the functions 

r* 

mo(t) = / p{t,x)dx, 

JXt 

ni-{t)- / p{t,x)dx, 

J— OO 
noo 

m+(t ) = I p{t,x ) dx, 

Jx* 

then mjitm) —► nij. 

Remark 5.7. Suppose that x_ = xq = x*, hence p m {mo + m-)S Xir + m + 8 x **. The lemma above proves that 
m + = lirn,„ m+{t m ), but does not gives any information about the masses m_/o. 

Anyway, since m- + mo + m + - 1 and for every t m the relation m_ [t m ) + rno(t m ) + m+ {t m )= 1 still holds, hence 
we can assert that 

m- + m 0 = lim [m_(f ra ) + m 0 (f m )], 

m 

which is going to suffice for the purposes of our analysis. ♦ 


After all these considerations we have in particular obtained that 

a{t m ) —■ H'lXi), 

which allows us to choose a triple (X_ ( t rn ) , Xq ( t m ) ,X + (t m )) m such that 

a{t m ) = H'{X i {t m )), (5.9) 

and Xi[t m ) —*• Xi by construction. Considering moreover equation 15.81 and Lemma fSlil we can hence say 
that 

Y m{tm)Xi{t m ) + 0{ 1), 

ze{-,0,+) 

with o(l) —► 0 as t m — ► oo. 

From now on we are going to refer as X, to functions that satisfy 15.91 , and not any more for solutions of 12.111 
with K — 0, +oo. 

Remark 5.8. a priori we could have different weak-* limit equilibria for p depending on the diverging se¬ 
quence of times. 

In other words, set {t m ), (t' m ) different sequences st t m , t' m —>■ +oo, set p m ,p' m the distributions, the above the¬ 
orem explains that p m — Z,p' m — Z', but might happen, a priori, that Z Z'. The next question is: is there 
uniqueness at the limit? If yes under which conditions? 


♦ 
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5.1 The uniqueness problem. 

In this section we want to identify some hypothesis under which the problem treated all along this paper has 
a unique weak-* limit as t —>• oo. 

The first step in our analysis is going to be the following lemma 
Lemma5.9. Consider D defined as in 15.41 . Than lim ,_ +no D(t) = 0. 

Proof. Thanks to Corollarv l5.5l we can argue that 

J [H\x]) 2 p(x, t) dx ^ max|iT '| 2 < oo. (5.10) 

With this consideration in hand, considering that D is defined as in 15.41 , it is easy to check that D, D' e 
Accordingly to 15.41 


D(f) = -| J H'{x)p{x, t) dxj + J [H\x)) 2 p(x, t) 


dxe L° 


= 0 - 2 (t)€l°° 


And 

D f [t) = -2o[t)- J H'{x)p t {x,t)Ax + J (H'(x)] 2 p(x, r)dx IbP+ ^^ 

2o{t)- J H"[x) [h\x) - cr(f)) p(x, t) dx + J (H'(x)) 2 p(x, t) dx ^ 

21|a|| L oo [msx\H'-H"\ + ||o1| L coj + max|H'| 2 < oo. 

Where in the last inequality we applied Corollarv l5.5l 

At this point, suppose it is not true that lim^oo D( t) = 0. Hence there exists a sequence such that D{t n ) > e for 
some £ > 0. Define the set A e = {t : D(t) ^ £}., and the set A f:jl — A e n [n,oo). Moreover A EiH = U kA E n , where 
A k E n are the connected components of A £i „, indexed by k. It is true that C —>■ 0 as n —► oo, otherwise D 

wouldn’t be L 1 (05, d t) contradicting 15.51 . this implies that also C j —► 0, for every index k. This considera¬ 
tion is valid for every rj > 0, in particular for el2 . 

Select a t n e A En ,s n t A E / 2 , n - These two sequences can be selected in such a way that \t n - s n \ - °°> 0. 
This is indeed true since (s n ) n can be chosen such that I t n - s „| sC 2 - diam ^ A k /2 n j - °°> 0, infering via La¬ 
grange theorem to state that there exist a sequence {T n ) n such that \D(t n ) - D(s n )\ - \D'(r n )\\t n - .s' n |, but since 
t n e A En , s,i € A £ / 2 , n is easily obtained that \D{t n )~ D[s n ) \ ^ |, but, considering that \ t n -s n \ - °°> 0 this would 
imply that | —► oo, contradicting D' e L°° and concluding the proof. □ 

This lemma states a very important property, which is that, for every diverging sequence of times we can 
extract a subsequence such that p m Z/e(- o,+} mu*> Xi ■ 

We would like to understand better the structure of the invertible branches^/ of H. 

Lemma 5.10. Assume the potential H is C* oc (D5), and consider a neighborhood [cr* where 8 is consid¬ 

ered to be small. Recall that, accordingly to 15.91 the functions X;, i - -, 0, + represent the invertible branches of 
the potential H. Then the functions A_/o (cr) = X-/o (cr) - x* are uniquely determined, A-/q (<j) = O [Vcr* -<jj in 
a vicinity ifcr*, and, moreover, setting a* -a = Act under the regularity assumption made on H, and supposing 
that H'" (x*) A 0 the following expansion holds 

A_ / 0 (cr)= Y, a~ /0 {Acr) n/2 + O[[Ao) 2 ). 

n =1 
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On the other hand, always in the same vicinity of ct* , denoting accordingly to the notation introduced, 
X+ (a*) = ct** we haveX + (ct*) - x** = 0(Act), and 

X + (ct*) - x** = £ a* (Act) n + O ((Act) 2 ) . 

72=1 

Moreover \ A-io (ct)| , |X+ (ct*) - x**| > 0 if a f cr*. 

Proof. The uniqueness is clear and comes from the definition in 15.91 considering that, by hypothesis, we have 
been considering potentials H with strictly monotone invertible branches. 

By definition of the function Xo (see 15.91 ) we have that a - H' (Xo (ct)), performing a Taylor expansion of the 
right hand side of this equation in terms of the perturbation Xo (ct) - x* and considering the fact that H" (x*) = 0 
we obtain that 

H'" f x 1 

cr — ct* =-(Xo(CT) - x*) 2 + O ((Xo (CT) - x*) 3 ). (5.11) 

Now, from equation 15.111 we can assert that H"'{x + ) f 0 comparing the signs of the left hand side with the 
right hand side, moreover, considering that by hypothesis H'" (x*) f 0 we obtain that H'" (x*) < 0. 

We need a detailed analysis of the factors A-/o, where A_/o(ct) = X_/o(ct) - x*. Let us make the following ansatz 

3 

A_ /0 (ct) = £ a- ,0 (Aay 1 ' 2 + O ((Act) 2 ) , (5.12) 

72=1 

which we will justify at the end of this proof, see Remark [5.11l Act = ct* - ct. 

From 15.111 , that 

A 0 (ct) = i /-(CT* - CT) 1/2 + o((CT* - CT) 1/2 ) , 

V |H"'(x*)| 1 

The same procedure gives us that 

A_(ct) = -\l ,„ 2 (ct* - CT) 1/2 + o((CT* - ct) 1/2 ), 

V \H' ’"(x*)r 1 ’ 

where indeed we have that A-[a) - X_(ct) - x*. What is left is to understand the asymptotic behavior of the 
linear term X+ (ct*) - x**, but this is easily obtained performing the same procedure above, in particular we 
obtain 

X + (ct) = X + (ct*) + -|^ + 0 (ct-ct*), (5.13) 

hence there is a linear dependence from the parameter ct, since we do know that H"(X + (ct*)) > 0. Equation 
15.131 can be justified formally via an argument similar to the one performed in Remark [5.11l 
Putting together the results obtained we get 

X+ (ct) = X+ (ct*) - h(u - ct*) + o (ct - ct*) , 

X 0 (ct) = x* + c(ct* -ct) 1/2 + o((ct* -ct) 1/2 ) , (5.14) 

X_ (CT) = X* - C(CT* - CT) 1/2 + o((CT* - CT) 1/2 ) . 


^ erec =\[w^T\’ b =W&o- 

At this point the first order expansion is clear. We will need though in the following the expansion of A_/o 
in 15.121 up to the linear term, i.e. the second order. 

To do so we evaluate the next Taylor element in 15.111 we obtain that 

Act= iH , "(x*)Ao(CT) 2 +iH t4) (x*)A 0 (CT) 3 + C)(Ao(CT) 4 ) = ciA 0 (ct) 2 + c 2 A 0 (ct) 3 + O(A 0 (ct) 4 ) , (5.15) 

indeed moreove we can express Ao(ct) as 

Aq(ct) = «J(ACT) 1/2 + fl2ACT + (9((ACT) 3/2 ), (5.16) 
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plugging 15.161 into 15.151 and after some algebraic manipulation we obtain 


Act = (a 3 ) 2 CiAcr+ C2 + ^2 c i) (A CT ) 3/2 + o((Act) 3/2 ), 

equating the coefficients of Act and (Act) 3/2 to zero and solving the non-linear system in the unknown we 
obtain two solutions, namely the two couples (aj, a®) and (rq , af ) 



We recall that c\, C 2 are defined as in 15.151 . 


(5.17) 

(5.18) 


□ 


Remark 5.11. We want to justify equation 15.121 . Taylor expantion yelds Act = cA_/o(ct) 2 + 0(A_/o(ct) 3 ) , with 
c > 0. For ct sufficiently close to ct* both right and left hand side of the equation above are positive, hence it 
makes sense to take the square root on both sides obtaining the following two equations 


(Act) 1/2 = A 0 Vc + O(A 0 ), 
(Act) 1/2 = -A-\Zc + <D(A-l 


Define F° ((Act) 1/2 , Ao) = Aq \J c + O (Ao) - (Act) 1/2 , a straightforward computation shows that ((Act) 1/2 ,0) = 
\[c / 0. By the Implicit function theorem we argue that there exists a function do = do ((Act) 1/2 ) such that, 
Act = cA 0 ((Act) 1/2 ) 2 + O (a 0 ((Act) 1/2 ) 3 ) , moreover Ao was a C* oc (R), hence we can express it as Ao ((Act) 1/2 ) = 

3 

a° n (Act)" /2 + O ((Act) 2 ] , which is exactly 15.121 . A similar approach is valid also for A_. ♦ 

n =1 


At this point we can study the problem of the uniqueness as t —► 00 for potentials H which are CjL (K). 

Proposition 5.12. Suppose that the potential H satisfies the hypotesis stated at the beginning of this paper, in 
Subsection \2.1\ and moreover suppose that alt) e (ct* + 7], a* - 77 ) for some rj > 0 and for all times t > to- Then 
the limit is unique. 

Proof. First of all we claim that, there exist a local maximum Molt) of p(x, t), around which all the mass of the 
unstable region concentrates such that, lim^oo | Mq ( t) - Xq ( t) | =0, where Xq ( t) is the only spinodal state such 
that alt) = H' (Ao(t)). Mo exists thanks to Theorem l5.2l 

Moreover we claim that Mq e Cf (R). This in indeed true since we are considering the mass transported 
along characteristics, and by Proposition l4.2l these are defined and C 1 globally in R. 

Since Mo is a local maximum for p satisfies p x (Mo(f), t) - 0, and hence, expanding the equation 1MSU into 
p t - H"lx)p+ ( H'lx ) - alt)) p x we obtain that 

Pt (Molt), t) = H"lM 0 lt))p (Molt), t) . (5.19) 

Set p IMolt), t) = rlt). Thanks to the hypothesis alt) e (ct* + q,a* - 77 ], and since Mq - Xq + o(l) we can state 
that Molt) e (x* + 77 ', x* - 77 ') for some 77 ' small and for t sufficiently large. By the structure of the potential, 
hence 

H"IMolt)) ^ -e,e>0, (5.20) 

hence, considering 15.201 we obtain r'lt) f -erlt). We can, at this point, apply Gronwall inequality to the 
previous inequality, obtaining rlt) ^ r(to)e _e ^ _to) 1 °°> 0, Proving that molt) —»■ 0 as t —► 00 . 

Moreover, under these assumptions, tfi+ > 0, hence 

e* = [m°° + 0 (D) X- (CT) + (m“ + 0 (D) X+ (ct) + o(1)X 0 (ct), (5.21) 
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(5.22) 


where m“ = lim^oo m+(t), but, by construction X,- e L°°, hence equation 1 15.211 can be restated as 

£* = m”X_((j) + m^X+ia) + o(l). 

We can, at this point, to prove the uniqueness. Suppose we do not have uniqueness for p, than there’s no 
uniqueness for cr either, which means that cr oscillates between two values cri < 02 , and hence at these two 
values the convex combination for £* in 15.221 would give two different results, in particular 

£* + o(l) = m°fX-(cr 1 ) + mffX+icr 1 ) < m°fX-(o 2 ) + m“X+(cr 2 ) = £* + o(l), 

which is indeed absurd being £* a conserved quantity. □ 

Remark 5.13. It might be interesting to notice that Proposition 15. 121 underlines clearly the equivalence be¬ 
tween the vanishing of the spinodal mass and the uniqueness of the weak-* limit. ♦ 

Theorem 5.14. Under assumptions of Subsection\2J\on the potential H, iff* e [x*,x*] c , ifH isC 4 aroundx* 
and x *, and (x*), (x*) > 0, the limit is unique. 

Remark 5.15. In the following proof o(l) is going to be a general perturbation depending only on t such that 

o(l) — 0. 

Proof. Define a = limsup f ^ 00 o’(f),o;= liminf t ^ 00 o’(f),and consider the case in which £* > x*. The other case 
is simply symmetric. Note that this in particular implies that a > cr* . 

We are going to divide the problem in several simple sub-cases. 

1. suppose a ^ a*. Than £* = X + (cr) + o(l), and hence, since £* has to be constant <7 has to converge to a 
unique limit. 

2. the case in which <7 < cr* has already been discussed in detail in Pronosition l5.121 

3. suppose er < cr* and 7f > a*. This implies that we can choose an arbitrary large time f* such that cr(f*) = 
cr* and such that there exist f* < h < t .2 such that a* <ti \ - (T(t\) < 02 - UU 2 ) < o 7 - 

Let be x\ = X + {(T\),X 2 = X + (cr 2 )>indeed if cri < a 2 by strictly monotonicity of H' we can say that x\ f. X 2 . 
We obtain hence 

£* — X\ + o(l), 

-x 2 + o(l), 

which is absurd. 

4. cr* < £ < cr* and a — cr*. In this case we are not allowed Gronwall inequality as in Proposition s. 121 but 
under these assumptions we know that m+ f 0, hence there exist 

lim m + (t ) = m°f. (5.23) 

t—*oo 

We consider some a close to the extremal value cr *, a < cr*, we are going to perform our analysis on the 
conserved quantity 

£* = m-{t)X- (cr(r)) + mo(t)X ,0 (cr(r)) + m < f’X + [a[t)') + o(l), (5.24) 

as long as cr / cr*. We will be forced to divide the proof in sub-cases. 

(a) Suppose fT"'(x*) < 0, whence in this case the asymptotic performed in Lemma S. 1 Ol holds. 
Inserting 15.141 into 15.241 . we obtain 

£* — mffX+lo*) + x* (m_(f) + mo(f)) + c(mo(f) - m_(f)) (Acr) 1/2 + O(Acr) + o(l). (5.25) 

We remark the fact that, thanks to 15.121 the term O (A cr) is Cj 2 c (K). 

Now, £* a constant it has to be independent from (Acr) 1/2 , to this end evaluate 

d£* 

° ~ d((Acr) 1 / 2 ) = c{m ^t)-m-U)) + 0({/\(7) m ). (5.26) 
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From the equation above whence we obtain that 


moU) = m_(f) + 0((Acr) 1/2 ). 


Moreover 

molt) + m.-(t) = 1 - m+[t) — 1 - + o(l) = B + o(l). 

Considering the result above with the one in equation 15.271 we can argue that 


ni-{t ) = 


£ + o(l) + e>((Ao-) 1/2 ) 


’ + o(l) + O ((Act) 1 ' 2 ) . 


Moreover a similar approximation is valid also for mo thanks to 15.271 . whence 

molt) = + o(l) + 0((Acr) 1/2 ). 


(5.27) 

(5.28) 


(5.29) 


(5.30) 


This approximation is valid as long as the approximation 15.251 holds. Thanks to 15.271 we can say 
that m 00 = m™, whence, if n oscillates, close to the bifurcation point the masses stabilize around 
the same value. 


Consider now the term c (mo (t) - m_ ( t)) appearing in 15.251 . We know, thanks to Remark 15, 1 61 that 
c(moU) - m-(t)) = -2 (aj /0 (1 - m”) - m™b) (Acr) 1/2 + £>((Acr)), (5.31) 

withe = |«7 /0 | > 0. Notice that thanks to 15.181 and the fact that H (4) (x*) > 0 results that -2 ( a 2 (] (l - m™) - nrf b) > 

0. 

At this point, consider an interval [ ti, ? 2 ] such that a < a* is decreasing in such interval. This in¬ 
terval always exists since a e C 1 and <r < a*, and such that cr (f 2 ) is close enough such that the 
approximation 15.311 still holds. Thanks to 15.311 we can argue that in [t \, t 2 \, mo - m_ is an in¬ 
creasing function. On the other hand as long as cr < a* we have that tho Sj 0 and /«_ 7 0, hence 
mo - m- has to be decreasing according to this consideration. This is indeed an absurd, proving 
that a can not oscillate. 

(b) Suppose at last that H'"(x*) = 0, whence, with the same considerations which have been done 
before 

* tf ( 4 ) (X*) 3 . 

a-a =--- A_ /0 (cr) + O (A 4 /0 (cr)) . 

Recall that we are considering H ( 4 ) (x*) > 0 . Set Act as in the point a. With a procedure similar to 
the one performed in Remark l5.11l we can conclude that A-i 0 {a) = c(Ao-) 1/3 + o((Acr) 1/3 ). 

( FT(4) fy \ _ 1/3 

c in the equation above that takes the following explicit expression c =- < 0 , whence, 

for cr < a* wehave that A_/o (cr)-x* = c(Act) 1 / 3 +o((A(j) 1/3 ) < 0 , This indeed implies that X_/o (cr) < 
x* for a sufficiently close to a* , but this contradicts the definition of Xq, hence we obtained an 
absurd. 


□ 


Remark 5.16. We want to justify equation 15.311 . Define R(cr) the 0((Acr)) function appearing in 15.251 . 
Whence 


d R 

d((Acr) 1/2 ) 


C>((Act) 1/2 ) 


where O ((Acr) 1/2 ] is the function appearing in 15.261 . Thanks to 15.121 and 15.131 we can express R{a), i.e. 


i?(cr) = molt)a°Aa + m-{t)a 2 Act - m^bkcr + o(l) + 0((Acr) 3/2 ]. 

Recall that, thanks to 15.181 a® = a 2 = a, and that the term C9((Arr) 3/2 J is C 1 around 0 and that o(l) depends 
only on t. Using 15.281 we obtain that 

R{cr) — (ct(l - m“) - m™b) Act + o(l) + 0((Acr) 3/2 ). 
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At this point differentiate both sides of the equation above obtaining 


which proves 15.3ft . 




- m'Jfb) (Ac) 1/2 + O (Act) , 


♦ 


A Estimates for the fixed point theorem and other technicalities. 

Lemma A. 1. Consider the operator C as defined in 13.221 . than C is invertible. 

Proof. Indeed equation 13.221 can be seen as 

C(pit) = (pit) - mKo*(pit) - (1- m)K+-k (fit), (A.l) 

for the following convolution kernels 

K 0 (t) = ae~ at x\ o,+oo)M> 

K+it) = be~ bt X{o,+cofit). 

We want to check the invertibility of C in 13.221 . 

To begin our analysis consider first the regularity of the left hand side in 13.211 . By definition we have that 
(p[t) - o[t) - do is C 3 , in fact inherits the same regularity of H around xq, in some interval of the form (-oo, 7'] 
for some T eU. Hence by the convolution structure of the equation JA.lt we have that the left hand side, and 
hence the right hand side of 13.211 is C 3 oc . 

Since JA.lt presents convolutions it seems reasonable to perform some change of variable in such a way to have 
our functions defined on the positive real line, hence apply the Laplace transform to obtain some information. 
Setting z = t-K for x > 0 equation 13.221 reads as 

J f»oo roo 

e~ aK (p)t-x) dx - i\ - m)b / e~ bK (pit-x)dx, 
o Jo 

which is again a convolution equation, and, moreover the convolution kernels didn’t change structure. We 
are performing an asymptotic analysis for t close to -oo, hence it seems reasonable, at least at the moment, 
to consider t bounded from above by some value to. Hence we can write t - to -x for x > 0, setting (pit ) = 
(pito-x) - y/ix), we can read the above equation as i jrix)-ma / 0 °° e~ ax y/ix-x)dx -(1 —m)b / 0 °° e~ bK y/{x-x)dx, 
i.e. Cy>ix) = y/ix) - mKo * y/ix) - (1 - m)K+ * y/ix), which is an equation in convolution form defined on the 
positive real line. 

Performing the substitution Wit) = 14/(to - x) — V(x), equation 13.211 turn into 

Cyrix) = Vix), (A.2) 

which has the same regularity of 13.211 but is defined on positive numbers. 

We can hence apply Laplace transform on both sides of 1A.21 obtaining Ly/iO) (1 - mLKo (0) - (1 - m)LK+ (0)) = 
LV (0), were Lyr, LKq , LK+ have respectively the domain 

DiLyr) = {0:Re0> a], 

DiLK 0 ) = {0 : Re0 > - a}, 

DiLyr) = {0:Re0> -b}. 


Hence we can express Lyr id) = LVi0)/C(0) as a meromorphic function defined on the half complex line D ( I.Ko) = 
{0 : Re0 > - a}. Our aim is to invert the term on the right hand side of this previous equation via inverse Laplace 
transform. 

In particular we want to show that we can express for this particular case the inverse Laplace transform as 
a residual evaluation, to do so first we have to prove that |Li//(0) < Ml\Q\ C for some c > 0, hence 


I Ly/m- 


If 

I Jo 


—f) r 

e y/ix)dx 


I/' 

|J0 


1 d 

- e 

6 dx 


-Ox 


y/ix)dx 


IbP 

< 


1 a 

--e~ 8x y/ix) 

V 0 


1 



I r°° 

/ e~ ex y/'ix)dx 
I Jo 
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Where in the last inequality we have proceed as follows. Consider -^-ipix) = -4^-</>(?) = it), and, since 
cr(f) = f H ! (x)d x R(x, t) d x with A’ defined in 13.31 we obtain 

^CT(f) = H"(X + (t)) [<j(t) - if'(X+)) - [H"(x+) + O(X+(0 - x + )]C(X + (t) - x + ) < c < oo. 


Now, we know that ip is continuous, we have to check were are localized the poles of Li// i.e. the zeroes of C to 
understand if we can indeed invert the operator C. 

Set the equation 


C(0) = 1- 


ma 
a + 0 


(1 - m)b 
b + 0 


After some algebra we obtain two roots 0i = 0, 02 = -(1 - m)a— mb < 0, which have real part strictly smaller 

than - a, in this way we obtain the following expression for i// 


i//(x) = 2n Y' Res 
t=l,2‘ 9 = B - 


(n x mey\ 

i C(0) )’ 


where this last equation is justified by fT] Theorem 16.39]. 


□ 


LemmaA.2. The equation CG(t) -8o(t-a), is solved for the function Git- a) where 

G(x) = [<5o Cx) + G r (x)] Xh+, (A.3) 

with G r (x) = [ci + c 2 e S2X ] j R+ (x), ci ^ O,0 2 < 0. 

Proof. At first I want to show that every solution of 

J C<p{t)- 0, for t< to, 

] lim tp[t ) = 0, ^’4) 


satisfies (pit) — 0 for t< to. 

Performing the usual substitution t- to -x we obtain that 1A.4I is equivalent to 

j Cy/ix) — 0, for x > 0, 

1 lim i//(x) = 0. 

[ x—►OO 

Applying Laplace transform we obtain that Li//(0)C(0) = 0 if I) - {Re0 > -a], but C is different form zero in D, 
which implies that Li// have to be identically zero in D. This means that Lipid) = / 0 °° e~ 8x ipix)dx — 0, in D, now 
we can always see Laplace transform of a function / as Fourier transform of an associated function, up to a 
constant i.e. 


LfiB) = JJ° e~ Bx fix)dx = J + °° e~ ilmdx [i H+ ix)e~ ReBx fix)] dx = 

this is a more flexible way to face the problem. Is easy to check that jr + i*)e~ ReB ‘ipi») - ipg e L 2 (K), hence as 
long as 0 e D we obtain that iFiipg) - 0 in L 2 . Fourier transform is an invertible operator in L 2 , hence ipo ~ 0 in 
L 2 , but since i p is continuous we obtain that ipg is identically zero and hence ipix) - 0 for x > 0. 

Consider now the equation 

C(Pit) = 8 a it). (A. 5) 

Since ip — 0 for t< a, with the substitutions 


t-a = x, 
( + a - z, 
(pit) = (pit-a) - 1 pix). 
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Equation IA.51 turns into the following 

J fX nX 

e ~ a (x-0(p(Qd(_ (1_ m)b I e~ b(x ~^(p(C)d( = do(x). (A.6) 

o Jo 

A straightforward application of transform methods may not be efficient, given the strong discontinuity 
presented in the problem given by the function do, in this spirit we try to substitute cp with a suitable decom¬ 
position that may lead to a problem sufficiently regular to apply the Laplace transform. 

To do so, consider <p as <p - do + G r . Inserting this function in such a form equation IA.61 reads as: 

£G r (x) - mae~ ax + (1 - m)be~ bx , 


where the member on the right hand side is suitable for application of transform methods. Applying Laplace 
transform to both sides we obtain the following equation defined in the half complex plane {Red > - a] 

ma (1 — m)b\ ma (1 - m)b 

eT^~~eTr) LGrm = eT^ + ~eTT’ 



which is equivalent, after some algebraic manipulation to the following 

6{ma+ (1- m)b) + ab B{9) 

1 0 1 2 + ((1 - m) a + mb)6 C(0)’ 


We have hence obtained that LG r can be expressed as a meromorphic function which has 2 simple poles 
located at Q\ - 0,02 < 0. Heaviside inversion theorem (see, for instance, (6] for a proof of such), can be applied, 
expressing G r as 


G r (x) 


Bid 1) c e, x Bid 2) c g,x 

C'(0i) C'(0 2 ) 


Ci + c 2 e d2X . 


We’ve obtained that if G(x) is solution to JA.51 we have that G(x) = [do (x) + G r (x)] Xr, with G r (x) = 0 as long 
as x < 0 and G r £ L°°(IR) since 0 2 < 0. 


□ 


In the proof of Proposition l3Hl we use a fixed point argument, i.e. we want to show that the operator T 
defined by mean of formula 13.251 is indeed a contraction between Banach spaces. To do so we require the 
following estimates. 

Proposition A.3. Let be W be defined by means of 13.231 . than we have that there exists some positive C = 
C{a, b, m,8) < oo uniformly in 8, and t< to sufficiently negative such that 

\W(t)\^Ce~ 2at (A. 7) 


Proof. Recall that W is defined as follows 

W(t) = m[-aU 0 [poKt,cr ) + p 0 (t,cr)) + (1 - m) [-bU+[p+){t,o) + p + (t,cr)) + 

+ maf {U 0 [poKt,a)-U 0 [pKt,K,a))Q'(.K)dK+m f [p{t,K, o) - p 0 (f,cr)) Q'[K)dK <323t 
Jo? Ju 

Recall as well that 

P0 = O(Y 2 ) 

p+ = o{y + ) 
p = e>|(^£r flr +y) 2 j 

and Yq, Y+, Y — o(l) as t —► -oo. 

Recall as well, as shown in Lemma 1231 that |Fo/+(t)| < e i-a+S)t ( } lence |p 0/+ (f)| ^ C//e 2(_f ' +A) r . A straightfor¬ 
ward computation shows moreover that \Uo[pot+]{t,o)\ f Cue 2( ~ a+fi] r . By the definition of p{K,t) we know 
that p{K,t) — 0[{X(K, t) - xo) 2 ) , and, since lim^ ±00 |XCRT, r)| = [X+(t)| < oo, we can say that the function 
X[-,f) is bounded as long as the extremal characteristics X+ exist. In particular since X+ — x± + Y+ with 

Y+ (t) < e Ga+S) t we can sa y 

1 

X+(t) ^ x + + -, 

1 

X-(t) ^ x_ —, 

2 
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hence we can bound uniformly \X{K, t)\ f x+-x_ + l, which implies that | -4= j P p{K, t)e~ K ~ d /C| sC C#e 2(_a+5) t . 
We point out the fact that Q'(fC) = -^=e _Jf2 . 

A straightforward computation as before shows hence that | -j= f u U o [p] (K, t)e~ K ~ d AT| ^ C//e 2(_n+5) t . At this 
point, considering W as in equation 13.231 we can hence assert that 

| W{t)\^Ce 2[ ~ a+5)t (A.9) 


for some positive C- C{a,b,m,8,H ) < oo uniformly in 8 and t< to sufficiently negative. 


□ 


Proposition A.4. Let W be defined via equation 13.231 , and let (p\,<p 2 e K M f Q (8). Consider H with the same 
properties as in Lemm dZfil We recall again that all the functions pi, Y), i — 0, + have a direct dependence on the 
parameter a, which is equivalent as having a dependence for a parameter <p e K M f 0 (8), than 

| W{t,<Pi) - W(t,<p 2 )| < e~ {2a+S) r -1|01 -• (A. 10 ) 

Proof The functions p/, Y), i — 0, + have a direct dependence onpe K M j 0 (8), in particular the dependence is 
given by 

Yoi + -Uoi+[(p\-Uoi + [poi + ], (A.11) 

y = Wo[0]-Wo[p], (A. 12) 

and the functions p,- are defined as in lA.til . 

We want to derive first some estimate for the element \Uq [(pi - p 2 ) (f) |, with Uq [/] (f) = f_ co e _a(f_s) /(5) ds, 
hence 


I r ^ 

\Uo[<pi-(p 2 ) (f)| = I e ~ a[t ~ z) [pi{z)-p 2 (z))d; 

| J—OO 

< K Mto (S)' e ~ at J_ 


t 1 

e -(«+5 ) z dz < - L 


a -(2a+d) t 


-{a+ 8) 


(5)’ (A- 13 ) 


this estimate will be useful in the following. 

Since H satisfies the same properties as the H considered in Lemma 1231 we can say that, for to sufficiently 
negative, po [t,(pi] = ft (YoU, </>,•)), where we decide to express ft ( To (t, (/>/)) in its integral form, i.e. 


R{Y 0 [t,(Pi)) = f 
Jo 


Y 0 a,(pi) h '"(xo + 0 


[Yolt,Pi)-0 2 df, 


since it gives a better idea of the regularity of the reminder. We can express the reminder in such a way since 
by our hypothesis in Proposition l3.il H'" is locally absolutely continuous. Since Yo is at least C 1 (see Lemma 
12.31 . then ft(Yo) is C 1 as well. 

Consider at this point the function ft(x) = / 0 A H (A ° + ^ (x - () 2 d(. R e C^ 1 , u Cj 1 , hence for x, y in a compact set 
K, |ft(x) - ft(y)| ^ L(/C)|x- y|. We claim that 


lim L(IQ - 0. 
£(/Q-0,0e/C 


(A. 14) 


where C is the Lebesgue measure on IR. 

Indeed, for locally differentiable functions g on compact sets 1C, we have that ||g|Lo,i ^ sup|g'(x)|, which 

K X€lC 

means we can bound the Lipschitz norm in term of the sup of the derivative in the compact set. This means 
that, if we prove that R' (x) —<• 0 as | x| —► 0 we prove IA.141 . 

A computation shows that R' (x) = x / 0 A H'" (xo + z) d z - Jq zH'" (xo + z) d z —— ^ 0, since H e Cf oc . 

Now, R (Yfit,(pi)) — po(f,(/),), this means we have obtained that 

\po[t,(Pi)-po[t,<p 2 ) \ ^Ut)-\Y 0 [t,(pi)- Y 0 [t,(p 2 )\, (A. 15) 


with L{t) = o(l) as t —► -oo. 
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Moreover, since Yq is defined via equation IA.11L we can infer that 


| Yo ~ Yq (f, 02) | ^ |Wo (0t-02) (?) | + |Wo (po (0i) ~ Po (^ 2 )) (f)|> 


obtaining, 

|po (f,0l) - Po (?,02) I < i(f)( I Wo (01 - 02 ) (?) I + I Wo (po (0l) - Po ( 02 )) (?) l). 

Hence, considering the term -aUo[po(<p)](t) + po(?,0) in 0.231 . with these estimates, we can argue that 

|-a (W 0 (p 0 (0i) - po (02)) (?)) + Po (f,0i) - Po (f,02) | 

< L(?) |W 0 (0i - 02) (?) | + (-< 2 + 1(f)) |Wo (p 0 (0i) - po (02)) (01, (A. 16) 

hence what’s left to understand is how to bound the term | Uq (po (0i) - po ( 02 )) (?) | • 

Repeating the procedure in equation I A. 1 51 we can argue that, there exists a 0 < q < 1 that we can make as 
small as we want since q — sup {L{t)}, such that, 

KU 

| Wo (po (0i j - Po ( 02 )) (?)| ^ t?|Wo (Fo(f,0i) - y o (f,02))| 

q |Wo (0i - 02) (?) - W^ (p 0 (0i) - Po ( 02 )) (?)| 

N 

< L 4" \K +1 (01 - 02) (?) | + q N \ur (po (0i) - PO (02)) (01 • 

n =1 


Is easy to check that 

/ 1 \ iV+1 

|W 0 W+1 (Po (0i) - Po (0 2 )) (?)| < ( 7 ^ 5 ) e2 '~ a+S) ' 

which in turn implies the following 

1 " 1C (» M - Po W) ml o 21 -* 8 ' — 0, 

if q < - a + 28. At this point we can hence say that 


|W 0 (po (0i) - po (0 2 )) (?)| < £ d"|W 0 " +1 (0i-0 2 ) (f)|, (A. 17) 

n =1 

which means that we have bounded the term \Uo (po (0i) - po ( 02 )) (?) | from above with another term which 
we can evaluate thanks to equation IA.13I . 

With the estimate in 1A.13I we can bound from above the term on the right hand side of 1A.17I with 

E r K*' (*r -<*) M IS E - p- 12 ""'' ■ 10. - 

which, considering IA.17I gives us 

l"° < w W w (<te » M| 4 ^ “' ,S) kr.« • (A ' 18) 

At this point we can hence plug the estimates in 1A. 1 81 and IA.13I into IA.16I . this gives us the following 


| - a (W 0 (p 0 (0i) - po ( 02 )) (?)) + po (?, 0i) - po (?, 02) | < 

ki {q, a, 8)e^ a+2S) f ||0 X - 0 2 1| k m tg (8) + fc [q, a, 8)e ~ {:2n+5) 1 1|0i - 0 2 II k m ^ (8) 

< k{q,a,8 , q)e~ {2a+5) r ||0i - 0 2 II jc m ^( 5 )• (A. 19) 

With the same procedure just performed, we can obtain the following bound 

|- b (W+ (p + (0i)-p+(0 2 )) (?)) + p+ (f,0i)-p+(f,02)| < k{a, 8 ,q)e~ ( 2 a+ 5 )t \\(pi-(p 2 \\K MtQ (S), (A.20) 
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and 



,-{2a+6)t 


1101-0211 K„ t m- (A.21) 


At this point, considering the bounds IA.19i - IA.21t and the structure of W which is given explicitly in 13.231 we 
can assert that 



\W(t,<pi) - W[t,(f>2) \ < e 


Concluding the estimate. 


□ 
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